Prom Inherent Structures Deformation to Elastic Heterogeneities 
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Using a well defined soft model glass in the framework of Molecular Dynamics simulations, the 
inherent structures are probed by means of a recently developed deformation protocol that aims to 
capture the Dynamical Heterogeneities (DH), as well as by the use of the isoconfigurational ensemble. 
Comparisons of both methods are performed by extracting the corresponding inherent characteristic 
length scales as the temperature of the system is cooled down from the liquid to the glassy state. 
The obtained lengths grow and depict an identical trend as the system falls out-off equilibrium, and 
appear to converge to the characteristic length scale that characterizes the Elastic Heterogeneities 
(EH) of the materials in the very low temperature limit, which is deeply related to the properties 
of the glass. This provides a first evidence of a relationship between DH and EH. 

PACS numbers: 61.20.Lc,64.70.P-,64.70.qj 



Nowadays, one of the challenging problem in the 
physics of glassy systems lies in bridging different con- 
cepts devised by recent theories of the glass transi- 
tion [ll [ij . One aim is to improve and test them, but also 
to develop new mathematical tools and frameworks. Re- 
cently, two fundamental concepts seem to be of predom- 
inant interest. They deal with addressing the "structure- 
to-mobility" relashionship by accounting for the involved 
collective phenomenae that appear as the glass transition 
is approached. 

To tackle this problem, most of the recent approaches 
develop a strategy in terms of monitoring a correlation 
length that has to grow as the glass forms. The na- 
ture, definition and pertinent choice of this length scale 
is still unclear From the "conventional" multi-point 
space-time correlation functions methods 043' ^ ""^y" 
namical" correlation lenght emerges that only accounts 
for the heterogeneous dynamics but fails to address the 
structure-to-mobility relashionship. On the other hand, 
several theoretical frameworks compose with a "static" 
length scale that should characterize t he g rowing emer- 
gence of spatially extent structures l3"13|, but its con- 
nection with cooperative dynamical events is still a mat- 
ter of debate. Furthermore, such a length scale is not 
confirmed experimentally and no consensus emerged re- 
garding its definition. Recently, one method was pro- 
posed IJ, |l5| , which aims to bridge the scale between in- 
herent structure formalism of supercooled liquids 1^ 121 
and the elastic properties of amorphous solids [l^ by ob- 
taining and monitoring a characteristic length scale. The 
obtained length was proposed to reveal the trace of an 
associated critical length diverging at the glass transi- 
tion, which in a theoretical framework derives from criti- 
cal phenomena and Random First Order Transition the- 
ories (RFOT) [l9l - [2^ . In the following, this formalism is 
applied to a recently studied model glass |23| for which 
"conventional" methods have been shown to not correctly 
account for the emergence of heterogeneities in the dy- 



namics and elastic properties of the resulting glass. 

Numerical model and Static length scale: The 
systems under considerations have been well charac- 
terized in previous studies (T^]. Briefly, they consist 
in slightly polydisperse two-dimensional Lennard- Jones 
glasses (2DLJ), composed of = 10000 particles at con- 
stant density, chosen in such a way that in the glassy 
state, the average pressure is close to zero, and that the 
glass mimics the behaviour of a soft glass. The criti- 
cal temperature at which the system falls out-off equi- 
librium has been evaluated around ksTc — 0.24e. The 
2DLJ glass, already in a ksT state X° = {r°}, 
is then submitted to a macroscopic elongation applied 



on each particles i, such that 
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rectly followed by a minimization of the potential energy, 
which gives rise to the final state r°''. Following the no- 
tation of Ref. [l3|, this defines the DQ deformation, here 
in the athermal limit. The noisy part of the displace- 
ment field (r°''^ — r"'''), i.e. non-affine component, has 
been shown to depict collective motions of particles or- 
ganized in large vortexes of characteristic size ^naf / 30 
interatomic distances [l3| that also determines the elas- 
tic heterogeneities (EH) length scale, and below which 
the classical continuum theory of elasticity is subject to 
strong limitations. 

Inherent structures via deformation formal- 
ism (IQD): A generalization of the above scheme to any 
initial ksT ^ state X = {r^} has been proposed in 
Ref. [31, which aims to capture dynamical "elastic" het- 
erogeneities. It consists in first generating the DQ defor- 
mation, finally written as X'^'^ = {r^^''}, while another 
path is performed by first minimizing the initial state 
X — {vi}, which gives the state rj, and then deform it, 
leading to the path X'''^ = {r^"'''^}. This second path 
is denoted as QD deformation. The comparison between 
the two DQ and QD paths d; = rf'^'^ —rj'^'^ finally defines 
the T-non-ajfine field. Both T-non-affine and non-ajjine 
displacement fields can be characterized by their mag- 
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Figure 1. (Color online) Magnitude of the T-non-affine at dif- 
ferent deformation strength e^x and temperatures fcsT. The 
field becomes independant of the strength from which it was 
generated at each temperature, below e^x = 2.5 x 10~^. 

nitude. In the ksT — >■ 0, the non-affine field has been 
shown to lie in the reversible deformation regime, as 
well as deformation rate independant for a range of val- 
ues 10^^ < e.j:x < 10^^. For the T-non-affine field, the 
measure of its magnitude d{kBT, e^x) = {N~^ d^)-^/^ 
as a function of the temperature ksT is a way to quantify 
the same rate effects. 

Results are summarized in the FiglU which shows the 
magnitude of the T-non-affine field at different deforma- 
tion strengths and temperatures. It becomes indepen- 
dant of the applied deformation, i.e weak perturbation 
regime, for all the considered temperatures for values be- 
low exx = 2.5 X 10~^. As in the case of the ksT 
non-affine field [2^, one can interpret this weak pertur- 
bation regime to the reversible path of inherent struc- 
ture deformations, namely that the energy landscape of 
local inherent structures is weakly deformed so that local 
structures can reversibly drop from one local minima to 
the closest neighbouring one, which is probably the same 
than the one before the perturbation is applied. For de- 
formation rates Cxx > 4.0 x 10~^, a minimum appears in 
the Fig [1] that also marks the onset of rigidity of the ma- 
terial when the temperature is decreased. Indeed, when 
the material starts to fall out-ofF equilibrium, one expects 
that it becomes more "rigid" and "responsive" to the ap- 
plied deformation. In the linear response regime, this 
should lead to a decrease of d{kBT,txx)- For stronger 
perturbations, more irreversible local paths arc explored, 
which lead to an additional "plasticity" and drives the 
"noisy" T-non-affine field to increase as fc^T decreases. 
In the following, we only consider exx < 2.5 x 10^^ and 
therefore writes the T-non-affine field like d{kBT). 

Inherent structures via isoconfigurational 
ENSEMBLE FORMALISM (ICE): The inherent structures 
are explored in the framework of the isoconfigurational 
ensemble (ICE) [I^l, which probes the dynamical "mo- 




Figure 2. (Color online) Spatial distribution of propensities at 
ksT — 0.42e in ICE formalism compared to the corresponding 
T-non-affine displacement field from IQD formalism. The 
same reference state X = {r;} is used. The colder the color, 
the lower the particle propensity, and vice versa for warmers. 

bile" heterogeneities of the material. Technically, we 
quench an initial state X = {vi} identical to the one 
used for the IQD, which leads to X"^ = {r'}. From this 
state, one performs Nr = 1000 runs over a simulation pe- 
riod of 1.5tq,(T}, where Tq,(T) is the a— relaxation time 
of the model [23|, while for each new runs momenta are 
randomly distributed from a MaxwcU-Boltzmann distri- 
bution. Then, the squared displacement (Arf)ic charac- 
terizes the propensity of motion for a particle i. 

Definitions: It is useful to separate the propen- 
sity domains in three parts: ft := 1 \ {{Arf)ic € 
[0.0, Ar|}. where Arf is the mean propensity magni- 
tude, n ■= 2 \ {{Ar^)ic e [Arf,w + Arf}, where w 
is the width of the distribution function of propensi- 
ties that gives the mean Ar|, and finally, := 3 | 
{{Arf) ic > w -\~ Ar^}. We therefore define the prop- 
erty Cn as the condition that must fulfill a particle i to 
be part of an ensemble f2 (bracketed part of the fi's def- 
initions), and we write the ensemble average for a quan- 
tity X, as {X)n = N^' E^^i HCn - X,) , while {X)n 
represents the normal average. Finally, in order to ex- 
tract local quantities from both IQD and ICE formalisms, 
one defines the one-to-one correspondance filter J-n := 
{{i,j) G n, AND, {di{kBT),dj{kBT)) e An} such that 
A n ■■= [{d{kBT))n - w g; {d{kBT))n + wn] where wn = 
y/var{P{d{kBT),n)}, and P{d{kBT),Q) is the one-to- 
one distribution function like the ones depicted in the 
insets of the Fig. [31 

Results and discussion: Fig. [2] compares the 
propensity map at fc^T = 0.42e > fc^Tc to the T-non- 
affine displacement field for the same initial configuration 
X = {vi] obtained using IQD. This highlights the cor- 
relations between the zones that follow the largest wave- 
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length excitation (negligible d{kBT) or weakly pertur- 
bated inherent structures) and the highly probable im- 
mobile paths (lowest propensities). It also shows that to 
the softer zones in the energy landscape (large propensi- 
ties) seem to correspond large fluctuating zones in terms 
of local inherent structure properties with larger d(fcBT). 
This is rationalized in the Fig. [3] where one plots the his- 
tograms of the one-to-one correspondance between the 
norm diikBT) and the corresponding (Ar^)ic for a par- 
ticle i. To get rid of the temperature dependance, the T- 
non-ajfine fields are rescaled according to their minimal 
and maximal amplitude, respectively Min^ and Max^, 
which are calculated over the whole set of N particles. 
As the temperature decreases, the distribution evolves 
from a flat one to an asymetric one, meaning that par- 
ticles with high propensities also have a large d{kBT). 
This is enhanced as the systems falls out-off equilibrium 
and confirmed in the insets of the Fig. [3J where the one- 
to-one distribution P {d{kBT), $7) is computed for parti- 
cles that belong to an emsemble f2. For the lowest tem- 
perature kBT = 0.42e, it clearly shows that particles of 
high propensities, i.e. := 3, have a broad distribution 
of d{kBT) that also strongly departs from the average 
{d{kBT)) N by larger values. The opposite trend is ob- 
served for particles belonging to the low-level ensemble 
VL := 1. The same distributions are also plotted at an 
higher temperature fc^T = 0.8e. One notes that they 
all look closer to the average {d{kBT)) m and that the 
one-to-one correspondance distributions per ensemble VL 
strongly overlap each others. At this temperature, the 
correspondances are mixed up and it is hard to correlate 
dynamically facilitated domains to the local structural 
and mechanical properties of the material. 

In the FigHI^a), the IQD formalism is used in order 
to obtain the "elastic" coefficients through the Hooke's 
law (T^]. This is done by estimating the difference in 
stresses, for the two DQ and QD branches, between the 
final X'^'' and X'^'^ and their respective initial states. The 
two elastic Lame coefficients derive from the relations 
= ^cT^y/e^, and = (Act^^ - Act^J /2e,„ where 
/3 = DQ or QD. The shear modulus /z^ {kBT) are plotted 
in the FiglUJa), for €xx = 10^'^ and for both /3 protocols. 
One notes that the QD branch is insensitive to the tem- 
perature, mainly because the inherent structures are first 
probed and then weakly deformed, which can only lead 
to an average shear modulus that reflects the averaged 
intrinsic material properties of the glassy state. This is 
confirmed when one compares (ksT) to the ather- 
mal limit Conversely, the DQ branch first deforms 
and explore the pre-probed inherent structures and the 
possible neighbouring paths, and thus reflects the emer- 
gence of the intrinsic properties of the glass as the tem- 
perature is decreased. fjP''^ {kBT) varies from negative 
(unstable) values to positive ones (resistance to deforma- 
tion) when fcsT decreases, and reaches (fesT) when 
/cbT « ksTc at which the system starts to fall out-off 
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Figure 3. (Color online) Main Panel: histogram of the one-to- 
one correspondances between the norms of the T-non-affine 
field for particles i and their propensity. Insets: one-to-one 
distribution function P {d{kBT),i}) of the T-non-affine field 
for particles i belonging to propensity ensembles Q, and for 
ksT > ksTc- The vertical dashed line corresponds to the 
average value {d{kBT))N (see text for details). 
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Figure 4. (Color online) (a): shear modulus obtained from 
the Hooke's law following the j3 = DQ or QD branches, as a 
function of the temperature and for a deformation rate exx = 
10~^. Its value in the athermal limit [l^l is also given, (b) and 
(c): obtention of the length scales and £,Hp^ that capture 
the spatial correlations in the T-non-affine field after use of 
the one-to-one correspondance filter at fcsT = 0.42e. The 
functions Cn and Hp^{Q.) axe described in the text. 



equilibrium. From this temperature, inherent structures 
are formed and selected, and define the local and average 
properties of the material. 

Making use of the seeming complementarity of 
ICE and IQD, we then employ the filter Jt2 in or- 
der to capture the distance r until which the fields 
d{kBT) are correlated. The use of the filter ensures 
that the one-to-one correspondance property is pre- 
served. One first considers the correlation function 
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Figure 5. (Color online) Temperature dependance of the char- 
acteristic length scales of the T-non-affine (a): from the 
exponential decay fit of B{h), and for two deformation rates 
in the weak perturbation regime, (b) : and €,Hp^ after use 
of the one-to-one correspondance filter Jri- For comparison is 
also given the ksT ^ non-ajfine characteristic length scale. 



Cn{Sr) = (d(fcBT, r + Sr) ■ d{kBT,r))Q which has been 
shown to correctly capture the structural dominant 

M, Sii. 



correlations in noisy displacement fields [18|, |2J] . Such 
a function is plotted in the Fig. IH^b) for each ensemble 
n, and from this are extracted the length above 
which the T-non-ajfine field is not correlated anymore. 
Alternatively, the vortex-like structures of d{kBT) 
and their location with specific propensity domains is 
quantified by computing the histograms of the second 
Legendrc polynomial function P2 = (1/2) ((cos^6') — 3) 
where 9 is the angle between two vectors, as a function of 
their distance from their origin, after applying the filter 
J"n- The characteristic length scale is therefore obtained 
by finding the distance r that obeys the property C-ffpj = 
{r I Hp^ (r, 17 1) n Hp^ (r, 17 2) n Hp^ (r, n := 3)}. 
The method is briefly illustrated in the Fig|4l[c) for 
the system at — 0.42e taken from the Figl^l 

Independantly to the one-to-one correspondance be- 
tween ICE and IQD formalisms, we also characterize 
the structure of the T-non-ajfine fleld using the 
coarse-grained correlation function [l^ [I^ defined as 



__,dj{kBT)eb{nj) 



B{b) = d{kBTr^{Eli [Nr'Ef^i-n--^J^»v^, 
where Otifij) = 1 if Vij < h, and otherwise, which 
is shown to be well approximated by an exponential 
decay of the form B{h) k, cxp[— and where £,b 
characterizes the extent of correlations in the field. 

Results for the temperature dependance of the derived 
length scales are compiled in Fig. [SJ The coarse-graining 
correlation length is also plotted for two deformation 
rates below t^x = 2.5 x 10~^. Its rate independence con- 
firms the weak perturbation hypothesis, and its tendency 
to reach the athcrmal value ^„a/ / also validates the abil- 
ity of the IQD formalism to capture the onset of coop- 



erative effects that emerge upon cooling and lead to the 
final structural properties of the glass. In the Fig. EJb) 
are plotted the length scales and S,Hp^ obtained af- 
ter use of the one-to-one correspondance filter J-q. They 
characterize the spatial extent of domains of propensities 
motion of particles, and how these domains depict a spe- 
cific inherent structure deformation ability, i.e. onset of 
mechanical response. A strong variation of both lengths 
therefore signals an increase of spatial heterogeneities in 
both an "elastic" and "mobility" sense. Such an increase 
upon cooling is indeed observed and, interestingly, the 
lengths appear to strongly increase around fcsT ^ feslc, 
and then slowly converge, for fc^T < /cbTc, to the ather- 
mal value £,naff that characterizes the elastic response 
of the glass. Hence, the lengths also mark the onset of 
rigidity of the material, to which corresponds a selected 
distribution of elastic heterogeneities that has emerged 
from the distribution of mobile ones upon cooling. 

Conclusion: In this work, the inherent structures 
of a model glass were probed upon cooling, using a re- 
cently proposed deformation protocol in conjunction to 
the so-called isoconfigurational ensemble. It has been 
shown that simultaneously using both methods provides 
a good framework that allows to estimate a characteris- 
tic length scale that can capture the dynamical hetero- 
geneities (DH), correlate them to the emerging mechani- 
cal properties of the amorphous material, and finally re- 
cover the elastic heterogeneities (EH) of the glass. 

However, the obtained length scales do not seem to 
behave like critical ones i-e. smooth temperature 
dependance, no divergence occurs, and no singularities 
appear around the Kauzmann temperature (extrapolated 
to fc^Tx — 0.162), such that one can not actually discuss 
the nature of the critical behaviour of the liquid/glass 
phase transition. Additional investigations are needed to 
clarify these fundamental points. 
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